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1. page 66, line 46 should read 
ENABC=Pl**2*PRESS*S0NIC*EPSL0*AP ARE/2 . /GAMA/SIG 

2. page 67, line 32, should read 
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FOREWORD 

The report presents the results of work accomplished at Pratt & Whitney 
Aircraft, Florida Research and Development Center under Phase III of Con- 
tract NAS8-21310, sponsored by the NASA Marshall Space Flight Center. All 
work was performed during the period 27 June 1970 through 29 February 1972. 

Mr. R. H. Counts was the NASA Technical Monitor and Dr. Uno Ingard of the 
Massachusetts Institute of Technology served as technical consultant. The fol- 
lowing personnel of Pratt & Whitney Aircraft contributed to the technical effort: 

P. L. Russell, J. W. Koenig, A. C. Schnell, and G. D. Garrison, who was the 
Program Manager. 
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Symbol Description 
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NOMENCLATURE (Continued) 

Description 

Effective Length Correction Factor, Defined by i~ t + h d 

Resonator Coefficient 

Absorbing Cross Section = W/l 

Frequency Ratio = f/f Q 

Array Specific Resistance Ratio 

Dimensionless Radiation Resistance = R /pc 
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Symbol Description 
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A. BACKGROUND ^ 

Past efforts at Pratt & Whitney Aircraft on absorbing liners have been 
limited to theoretical analysis and experimental evaluation of designs consisting 
of arrays of Helmholtz-type resonators. These devices are efficient combustion 
instability suppressors primarily because, when properly designed, they can 
absorb a large percentage, if not all, of the pressure wave energy acting on the 
surface of the liner. An array covering a significant portion of the combustion 
chamber surface is desirable to ensure stable combustion. However, in many 
instances (particularly cases of existing thrust chambers in which attempts to 
uprate performance are thwarted by the appearance of combustion instability), 
surface areas are limited and only a few resonators can be installed. Under this 
circumstance, the resonator array theory is not valid; in fact, both the absorption 
coefficient and the open area ratio used in the theory can only be arbitrarily defined. 

B. OBJECTIVE AND APPROACH 

For acoustic energy absorbers other than arrays of resonators, such as 
single Helmholtz resonators, slots, and quarter-wave tubes, a design analysis 
based on total energy absorption was believed to be a valid approach, but no 
theory was available for use in the nonlinear acoustic regime * where rocket 
chamber acoustic devices must operate. Effort under Phase III of the program 
was therefore devoted to the development and verification, using cold-flow acoustic 
devices, of such a design theory for nonarray absorbing devices. Other cold-flow 
acoustic experiments were conducted to investigate the dimensional limitations of 
acoustic suppression devices, including resonators, slots, and quarter-wave tubes. 
The acoustic data were obtained using existing impedance measuring devices at 
Pratt & Whitney Aircraft and from a new acoustic reverberation system that was 
procured specifically for use under the program. 

C. CONCLUSIONS AND RECOMMENDATIONS 

It is concluded that individual acoustical devices can be designed for maxi- 
mum acoustic energy absorption using the theory developed herein. Furthermore, 
the theory may be used to design multiple elements with arbitrary spacing between 
the apertures, but if the spacing is closer than a certain limiting value, the 
elements must be treated as an array. However, at this limit, no discontinuity 
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in the theory exists; the nonarray dependent variable, the resonator coefficient, 
is equal to the array dependent variable - the absorption coefficient. 

For the design of absorbing devices the following values of the resistance 
coefficient, K, are recommended: 

Resonators 0.37 

Quarter-wave tubes 0.75 

Quarter-wave slots 0.53 

Whenever possible, it is recommended that resonators be used in lieu of 
quarter-wave length slots or tubes because resonators inherently dissipate more 
acoustic energy per cycle, and they can be designed for wide-band absorption to 
allow for inaccuracies in predicting the frequencies of resonance and the fre- 
quencies of the oscillations to be damped. 

Quarter-wave length absorbing devices have fundamental frequencies that 
correspond to lengths greater than the actual. The effective length and the actual 
fundamental frequency can be approximated using a simple technique discussed 
herein. 

v The number of apertures used to produce a resonator of given total open 
area with a common, unpartitioned cavity is immaterial unless the device is 
cavity- volume-limited. If so, highest absorption will be obtained with the 
fewest number of apertures . 

The apertures of absorbing devices may be placed adjacent to or in the 
corners formed by injector faces, baffles, etc. , .with no detriment to the acoustical 
characteristics of the device. In fact, the power absorption will be slightly 
improved in a predictable manner. 

It is recommended that single resonators be designed so that the ratio of 
the aperture diameter to the product of the quarter-wave length and cavity backing 
depth is less than one. Otherwise, the absorbing characteristics of the actual 
device can be less than theory predicts. 

From experiments with resonators having long apertures with respect to 
the quarter-wave length, it was shown that absorption at the harmonic frequencies 
can be as high as that of the fundamental. An analytical technique for predicting 
the resonance frequencies and the absorption at any frequency was developed and 
included in the nonarray device design computer program. The program, which 
can be used to design both resonators and acoustic cavities, is included herein 
as Appendix C. 
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SECTION II 

DERIVATION OF THEORY 

For the design of nonarray absorbing devices for use in rocket thrust 
chambers, a comprehensive theory is required to relate the acoustical character- 
istics to the physical dimensions of the device. Other factors that must be in- 
cluded in the theory are the properties of the media in which the device operates 
and the effects of media flow, if present, on the acoustics of the system. A 
great number of publications dealing with nonarray acoustic elements for archi- 
tectural applications (especially single Helmholtz resonators) have been published. 
However, none has presented a theory suitable for rocket application, primarily 
because the sound pressure amplitudes and flow fields peculiar to rocket motors 
are rather unique. 

A. BASIC DERIVATIONS 

For a resonant acoustic device, e. g. , a single Helmholtz-type resonator, 
with dimensions small with respect to wavelength, mounted in an infinite wall 
subjected to sound waves of normal incidence, the total specific impedance is, 
by definition, the following function of the complex ratio of incident pressure, 


Pp to particle velocity, u: 

Z = 2 P j/u (1) 

If p and c are the media density and sonic velocity, respectively, the incident 
energy intensity is 

I = Pi 2 /2pc ' . (2) 

and the power absorbed by the device is 

W = u 2 RA/2 (3) 

where R is the dissipation resistance and A is the aperture area. The absorbing 
cross section, e', is, by definition (Reference 1) 

€' = W/I (4) 

which, by using the preceding equations, can be expressed as 

«’ = 4pcAR|z| -2 (5) 


As shown in Reference 1 the maximum absorbing cross section for a particular 
frequency of wavelength X Is 

«'max =x2 / 27r 
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If a resonator coefficient is defined as e = t 'A ' max , equations (5) and (6) may 
be used to produce 

e = 8irpc AR (a |Z| )“ 2 (7) 

For a single resonator in an infinite wall the radiation resistance. Reference 1, 


R r = 27rpcAX~ (8) 

which, when substituted in equation (7) gives 

C = 4 R f R | Z | — 2 (9) 

The significance of t, the resonator coefficient, is similar to that of the 
absorption coefficient used for resonator arrays. It is both a measure of the 
energy absorbing characteristics and the efficiency of a resonator. The resonator 
coefficient can only vary between 0 and unity; values near zero mean that essen- 
tially all of the acoustic energy incident in the vicinity of the resonator is being 
reflected (or reradiated) and little absorbed. Values of € near unity mean that 
the resonator is absorbing as much of the incident energy as can be absorbed. 
Values of e less than unity indicate that at least theoretically it is possible to 
select a different resonator configuration that will absorb more acoustic energy. 

The total specific impedance of an individual resonator is 

Z = R + R r - iX (10) 

which can be used to put equation (9) in the form 


9 9 vxx; 

(0 + ly + (X/R r ) 

where 0 = R/R r . From previous work, Reference 2, it is known that the dissipa- 
tion resistance can be computed using 

R = KpuC f -2 (12) 

where K is the resistance coefficient (0.37 for circular apertures), and Cj is 
the aperture discharge coefficient. 

For a Helmholtz -type resonator the reactance for use in equation (11) 
is found from 

X = -(A/A c )pc cot (wL/c) + wpi (13) 
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where u) is the angular frequency, A c is the cavity cross-sectional area, L is 
the cavity backing depth, and i is the effective aperture length. For high sound 
levels a good approximation for l is known from Reference 1 to be 

& — t + 6d 

where 

5= 0.85 (1 - 0.71 JT e ). (14) 

In the above equations t is the aperture length, d is the aperture diameter, and 

a , the cavity open area ratio is defined as 
c 

<S = A/A . 
c c 

B. MULTIPLE ELEMENTS 

The nonarray acoustic theory will now be extended for application to 
multiple elements with arbitrary spacing between the apertures. It will be shown 
that if the spacing is closer than a certain limiting value (which is a function of 
frequency) the elements can be treated as an array. Furthermore, at this limit 
no discontinuity in the theory exists; the nonarray dependent variable, the res- 
onator coefficient, is equal to the array dependent variable, the absorption 
coefficient. 

Let the absorbing device consist of an array of resonators with circular 

apertures of diameter d spaced a distance b apart in a square pattern. Although 

2 

the facing (or incident) area per resonator is always b , the maximum absorbing- 
2 

cross section, A /27T, can only be obtained if the facing area for each resonator 

2 

in the array is a least equal to A /2 t ; this means the limit in the single resonator 
theory occurs when 

b 2 = A 2 /27T (15) 

or 

b~0.4A (16) 

If b< 0.4A the maximum absorbing cross section is simply b 2 , and the resonator 
coefficient by definition becomes 

€ = e'b -2 = W/Ib 2 - (17) 
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where W is still as defined by equation (3). The quantity W/b 2 in the above 
equation is actually the absorbed energy intensity, I a , which causes equation (17) 
to become by definition the absorption coefficient, a of the array; i.e. , 


, _ w 

€ = — T 

Jb Z 


I 


= a 


(18) 


For the special case b = 0.4A, substituting equation (15) in equation (8) gives 
the expression for the radiation resistance for arrays, 

R r = peer (19) 

2 

where a is the facing open area ratio, A/b . Equation (19) is also correct for 
all values of b < 0.4A. 

To determine the absorption coefficient, equations (12) and (13) are used 
for R and X and the coefficient is computed from 


40 


(0 + l) 2 + (X/apc) 2 


(20) 


where 0, the specific resistance ratio, is defined as R/pca. 


Note that for b = 0.4A either equation (11) or (20) can be used to compute 
the dependent variable because 


R r = 2 7rp cAX -2 = pea 


( 21 ) 


which means that for this case € = a, and that no discontinuity in the theory exists. 
C. OPTIMUM CONDITIONS 

It is known for arrays that the condition a = l maximizes the power 
absorbed. Proof that a similar condition exists for individual elements follows: 


Using equation (12), the power dissipated, equation (3) may be written as 

( 22 ) 


W = l/2Apc 3 B" 2 m 3 


where B and m are defined as 


B = KC 


-2 


and 


f 


m = Bu/c 


(23) 


(24) 
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The relationship between particle velocity, incident sound pressure, and 
impedance at resonance is 

u 0 = 2 Pi /(R r + R) 0 (25) 

which, by using the above definitions can be written as . 

-1 

(26) 

If P is the static pressure, the perfect gas equations for density and sonic velocity 
can be used to obtain 
2 



pc = 7P 

where 7 is the ratio of specific heats. 


(27) 


By letting 0^ = R^/pc, using equation (27) and rearranging, equation (7) 


becomes 


2 2p i B 

m o + *r m o = IP 


(28) 


which has for a solution 


0 r 

m o = -— + 




2p.B 
P" 


It may be shown that for this case, 

0 =k 2 A/2 ir 
r 

thus, by substituting equations (29) and (30), equation (22) becomes 

, . 3 

W q = l/2pc 3 B _2 A 



i 2 a\ 2 2p i B k 2 A 

LV \ 

4 717 IP 47 T 


(29) 


(30) 


(31) 


This equation gives the power dissipated at resonance by a nonarray absorber as 
a function of aperture area and incident sound pressure. For design purposes a 
more useful expression would be one from which the optimum aperture area could 
be determined. Such an expression is obtained by equating with zero the first 
derivative of equation (31) with respect to A. The resulting equation, 


opt 


27T 
, 2 


7 P 


(32) 
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which may also be written 



(33) 


shows that for a given media the aperture area that produces maximum power 
dissipation is a function of both the frequency and amplitude of the incident wave. 


Substituting the above equation back into equation (31) gives the expression 
for the maximum power that can be absorbed, i. e., 


W 

max 


(A p i) 2 


(34) 


Using equations (4) and (6) with the definition of the resonator coefficient, 

e = e'/i 1 , produces 

max ^ 



(35) 


which shows that for a given sound pressure and frequency 6 will be maximized 

when W = W . The proof is completed by substituting equations (2) and (34) 
nicix 

into equation (35) which results in ^ 


= 1.0 


w 


max 


(36) 


Inspection of equation (11) shows that 6=1 can only result if X = O and 0 = 1 

i. e. , if R = R at resonance, 
r 

D. OFF-RESONANCE THEORY 


To obtain an expression for the off-resonance behavior of the resonator 
coefficient the quality factor Q, which is defined (Reference 4) as 

Q = W Q M/(R+R r ) 0 (37) 

is used; this parameter is a measure of the sharpness of resonance (or band- 
width absorbing characteristics) of a driven resonator. For a Helmholtz-type 
resonator the inertance is 

M = pi (38) 

thus, equation (37) can be written as 

Q= W o p|/(R+R r ) o (39) 
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or, using 0= R/R r and equation (8), as 


Q = 


* o l 


M0 O + l) 


(40) 


It is reasonable to assume that at frequencies near resonance the resistance is 
constant; then the total impedance can be written as 


Z = (R+R r ) 0 [l-iQ(f-l/f)] 


(41) 


where { is the frequency ratio, f/f . 


Substituting equation (41) into equation (9) produces an expression for the resonator 
coefficient as a function of frequency 


4/3 


e = 


(0 o+ i) l 2 i + q 2 (r-i/n 2 


(42) 


Where required in any of the above equations, the frequency of resonance, f , 
is obtained by solving equation (13) for X = O, or, if L <<A computed from 

1/2 


f o 2 7 r 


Vi 


(43) 


where V is the volume of the resonator cavity. 
E . PAST- FLOW THEORY 


Rocket chamber acoustic devices must usually operate with a net flow of 
combustion gases past their apertures . It has long been known that the past- 
flow causes the acoustic resistance of the absorber to be greater than if there 
were no flow and causes resonance to occur at a slightly different frequency. 
Attempts to quantitatively determine the flow effects have been made (Ref- 
erences 2, 5, and 6); however, results from the different experiments were 
inconsistent and, in addition, because of the limited range of velocities, had 
to be extrapolated for most rocket chamber applications. Recently, (Ref- 
erence 7) both cold flow experiments and firings of a dynamically unstable 
rocket motor containing small absorbers were conducted and the flow effects 
were measured with Mach numbers up to 0. 73. It was found that 

1. Effective aperture length with flow past the apertures is best described 
by 

l = t + 0.375 [0. 85d (1-0. lyfo)] (44) 

P c 
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where <7 is the cavity open area ratio, A/A c . This value 
of ip should be used instead of that from equation (14) to 
compute the frequency of resonance, reactance, and other 
design parameters. 

2. Acoustic resistance, if 



where Mp is the Mach number of the past-flow and T c and T a 
are the chamber and aperture gas temperatures (°R), respectively, 
should be computed from 


R p - T Bpc M p <T,/T a ) 2 + i- Bpu (46) 

If the above inequality (45) is not satisfied the effects of flow 
can be ignored. 


An expression for the optimum area with high velocity past-flow can be 
obtained by using the required conditions for t = 1, i. e. , R^ = R^, and X = O. 
For this case 


% (R + R ) 

' p r'o 


Setting equation (46) equal to equation (8), substituting equation (47), and re- 
arranging produces 

\ _ 2 7 r 

optjp k 2 
o 

For the special case 



either equation (12) or (46) can be used to determine resistance and either 
equation (32) or (48) for the optimum area. 



'BMp/ 6) 2 (T c /T a ) 4 
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F. THROUGH-FLOW THEORY 


_If there is net flow of Mach number Mp through the absorber aperture 
such that M t > u/2c, then the resistance is (Reference 2) 


P 



(50) 


and the optimum area becomes 



(51) 


For high incident pressures (160 db* and greater) the effects of through-flow on 
reactance (and therefore on effective length and resonance) were found to be 
negligible (Reference 2). Furthermore, if < u/2 the effect of through flow 
on resistance can be ignored. 

In some applications there will be net flows both through and past the 
absorber aperture; however, as shown in Reference 8, the absorber can be 

treated as if only through-flow were present. 

\ 

G. LIMITING CONDITION 

All of the above analysis is valid only if the viscous resistance of the 
absorbing device is small in comparison to the nonlinear resistance. From 
previous experiments with circular apertures (Reference 2) it is known that 
this condition is fulfilled when 

Pj/P > 6667 (c C f ) -2 (52) 

Rocket chamber acoustic devices are typically designed for dynamic pressure 
amplitudes of at least 1% of the static combustion pressure, i. e. , p^/P = 0. 01, 
which is well above the limiting value of equation (52). 

H. NONARRAY DEVICE COMPUTER PROGRAM 

A FORTRAN IV design computer program is included herein as Appendix C. 
For a device of given dimensions the impedance, resonator coefficient, and other 
acoustical parameters are computed at resonance and, if desired, over a range of 


*In this report all sound pressures expressed in decibels (db) are based 
on a reference pressure of 0. 0002 microbar. 
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frequencies and incident-to-static pressure ratios. The theory used in the pro- 
gram is identical to the above, with one exception: it was not necessary to assume 
the resistance was constant for the off- resonance computations because an iterative 
solution of the impedance equation (1) was used to determine the particle velocity, 
i.e. , u is computed by solving 
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SECTION III 

EXPERIMENTAL PROGRAM 
A. REVERBERATION SYSTEM 

A reverberation system was fabricated and used to provide cold-flow 
acoustic data for verification and improvement of the theory for the design of 
nonarray absorbing elements. The system as illustrated in figure 1 consists of 
a chamber, embedded inside a gravel-filled bunker, driven by two electropneumatic 
transducers mounted on exponential horns. The internal volume of the reverber- 
ation system, including the horns, is approximately 35 ft 3 . The chamber 
(figure 2) is pentagonal in cross-section with walls 2 ft wide and an overall 
chamber length of 4-1/2 ft; one endwall is tilted 10 deg off- vertical to prevent 
the formation of a uniform longitudinal mode. Gaseous nitrogen used to drive 
the electropneumatic transducers exhausts through four 1-in. pipes on which 
low-pass filters are mounted to reduce the amount of sound energy escaping 
from the chamber. All instrumentation leads are introduced to the inside of the 
chamber through a 2-in. pipe mounted on the sidewall of the chamber. This 
pipe was plugged at the chamber end so that it would not affect the results of 
reverberation tests. A l-l/2-in. pipe was provided to introduce a secondary 
gas flow for tests requiring a flow of gas past or through the acoustic element 
being tested. Large plates mounted on the side and endwalls of the chamber serve 
as compactors to help settle the gravel tightly around the chamber walls, thus 
reducing vibration and loss of sound energy from inside the chamber. 

Electronic equipment used in the reverberation measurements is shown in 
the block diagram presented as figure 3. A General Radio random-noise generator 
was used to provide a noise signal with a Gaussian amplitude distribution, and a 
Technical Products wave analyzer was used to select a narrow band of the noise 
(250 Hz) for use as a signal to drive the reverberation chamber. The narrow- 
band noise signal was then amplified and introduced into the reverberation chamber 
through two Ling electropneumatic transducers. Four Bruel and Kjaer l/4-in. 
microphones were suspended inside the chamber with the outputs connected to 
a four-way selector box. Each microphone signal could then be individually 
selected and passed through a General Radio wave analyzer; the Altered signal 
was recorded on a Bruel and Kjaer graphic recorder in the form of decibels 
versus time. Two microphone signals were connected to the phase-measuring 
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system so that pressure-phase data could also be obtained. A discussion of 
reverberation techniques that were used with the apparatus appears in 
Appendix A. 

1. Preliminary Experiments 

Cold flow tests were conducted to evaluate the system response, i.e. , 
reverberation decay times, with different types of test signals used to generate 
the sound field within the chamber. A random noise signal with a uniform 
spectrum density (white noise) iri the frequency range 20 Hz to 20 kHz and 
having a Gaussian magnitude distribution was used in the first set of tests. 



Figure 1. Reverberation System FD 44696 
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FE 100996 


Figure 2. Reverberation Chamber FD 44697 

Maximum sound pressure levels (SPL) of less than 135 db were measured at 
selected frequencies in the range of 500 Hz to 2000 Hz. By turning off the 
random noise signal it was found that the background noise resulting from the 
gaseous nitrogen flow through the pneumatic sound drivers was fluctuating with 
SPL's between 100 and 110 db. The difference between the test signal and the 
background noise in each test was less than 25 db, which resulted in short 
fluctuating decay curves from which accurate decay times were difficult to 
measure. 

A pure tone sinusoidal test signal was used in a second set of tests in an 
attempt to produce higher sound pressure inside the chamber. It was found that 
sound pressures greater than 170 db could be generated at frequencies up to 
1000 Hz with levels greater than 150 db up to 2000 Hz. With pressure amplitudes 
of this magnitude it was possible to record decay rates over a 30 db range completely 
within the nonlinear sound pressure regime. The disadvantage of using a pure 
tone signal is that organized modes were established inside the chamber resulting 
in a spatially nonuniform sound field with sound pressure variations of up to 
±12 db. 
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Figure 3. Block Diagram of Electronic Equipment FD 61300 
Used in Reverberation System 

Use of the Sabine reverberation method to. measure the absorption of an 
acoustic specimen requires that the distribution of sound energy be uniform 
inside the chamber. To improve the distribution of sound energy, a modulated 
frequency (warble tone) test signal was used in a series of additional cold flow 
tests. The modulated frequency signal was generated by using a Wavetek 
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sweep/trigger Voltage Controlled Generator (VCG). The modulated frequency 
signal is generated by selecting a mid-frequency value for the primary generator 
output and simultaneously applying a second modulation frequency parallel to the 
primary generator output. The amplitude of the modulation frequency signal 
controls the upper and lower frequency limits between which the primary generator 
oscillates. 

Analysis of tests made with the frequency modulated signal showed that 
a significant improvement in sound distribution was achieved; the fluctuations 
were reduced to a level of ±2 db. It was decided that this energy distribution 
was satisfactory, and that cold flow tests could be conducted to evaluate the 
accuracy of the reverberation system by measuring the absorption coefficient 
of liner arrays having known acoustical properties. 

2. Evaluation Experiments 

The evaluation tests were conducted using a modulated signal with a mid- 
frequency of 1000 Hz and a modulation frequency of 10 Hz with a bandwidth 
having the same peak width as the one-third octave filter (figure 4) used to 
condition the output signal from the microphones mounted inside the chamber. 

In each of three series of tests reverberation decay curves were obtained using 
twelve different microphone positions inside the chamber. 



Figure 4. Comparison of Test Signal Power Spectrum DF 90636 
and 1/3 Octave Filter Bandwidth 
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The first experiments were made without a liner inside the chamber to 
determine the average absorption of the chamber walls and to obtain the baseline 
decay rate; a total of 150 decay curves were recorded. In the second series of 
tests a liner with known absorption characteristics was mounted flat on the 
sidewall of the chamber. The liner was a full array having an incident surface 
area of 144 sq in. with an open area ratio of 5. 5%, aperture diameter of 0. 093 in. , 
liner faceplate thickness of 0. 1 in. , and a backing cavity depth of 1.5 in. 
Theoretically, the liner was tuned for resonance at 1000 Hz, which would produce 
approximately 1. 0 Sabine of absorption within the modulation frequency band- 
width. (See figure 5. ) A third series of experiments was conducted to determine 
if small absorption values (approximately 0. 5 Sabine) could be measured. In 
these tests the same liner faceplate was used but the cavity backing distance 
was changed to 0. 88 in. increasing the resonant frequency to 1350 Hz; as shown 
in figure 5 the average theoretical absorption in the 1000 Hz modulation frequency 
bandwidth was reduced to approximately 0. 6 Sabine. A total of 132 decay curves 
were recorded. 



Figure 5. Theoretical Absorption Characteristics DF 90637 
of Liners Used in Evaluation Test Series 

Analysis of the data was based on the proposal of Sabine that the absorption 

coefficient of an acoustic specimen be calculated from the equation 


0.92V 


a 


Sc 


(D'-D) 


18 


(54) 




Pratt & Whitney Rircraft 

PWA FR-4993 


where: 

V c = the volume of the reverberation chamber 
S = the surface area of the test section 
D' = the decay rate and 
D = the baseline decay rate. 

The decay times for a 20 db drop were measured from the reverberation decay 
curves for all three series of tests and then used to determine the decay rates. 
Typical decay curves are shown in figure 6 and a summary of the test results 
is presented in table I. 


Table I. Reverberation Test Results 


Type of 
Test 

Number of 
Data Points 

Avg 20-db 
Decay Time, 
sec 

Avg Measured 
Absorption, 

% 

Theoretical 

Absorption, 

% 

Experimental 

Error, 

% 

Baseline 

150 

0.28 

- 

- 

- 

High Absorption 

120 

0.20 

87 

99 

12 

Low Absorption 

132 

0. 23 

51 

62 

18 


For the resonant (high absorption) tests the agreement between theory and 
experiment is probably as good as can be expected for this type of experiment. 

The more significant difference in absorption between theory and experiment in 
the off-resonance (low absorption) tests is believed to be the result of the rapid 
change in absorption with frequency over the range of input frequencies (figure 5), 
thereby permitting a small error in the measured decay rate to produce a large 
uncertainty in absorption. From the above results it was concluded that the re- 
verberation method is the most accurate when the liner elements are tested at 
resonance, giving maximum absorption and a minimum change of absorption with 
frequency over the band of frequencies covered by the modulated test signal. 

/ 

However, the acoustic devices to be used to develop the nonarray theory would have 
have to be tested not only at resonance, but over a wide range of frequencies. 
Because use of the reverberation method alone would introduce considerable 
experimental uncertainty, it was decided to instrument all test sections for 
pressure-phase measurements and use the equations of Appendix B for deter- 
mining impedance. The reverberation chamber, driven with a sinusoidal 
signal would be used as a source of high-amplitude incident sound waves. 
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B. VERIFICATION EXPERIMENTS 

Acoustic experiments were performed to obtain data for verification and/or 
improvement of the nonarray design theory. The reverberation chamber was 
used as a sound source and the pressure-phase technique was used to obtain the 
impedance data. The nonarray devices used in the experiments were single 
aperture resonators (figure 7) fabricated from sections of pipe and a steel plate. 
The backing cavity depth could be adjusted from 0 to 7 in. by moving the piston- 
type rear wall"; an O-ring was used to acoustically seal the resonator cavity 
volume and to hold the> piston securely in position while testing. The resonator 
was located in the center of 1-foot square plate to simulate a resonator mounted 
in a wall. 

Four series of resonator tests were conducted. In the first series of 
tests the resonator had an aperture diameter of 0. 75 in. and a backing cavity 
depth of 0.78 in. producing a 750 Hz theoretical frequency of resonance. The 
resonator was placed on the bottom of the reverberation chamber and a pure 
sinusoidal signal was used to generate the sound field inside the chamber. The 
actual frequency of resonance was determined by varying the frequency of the test 
signal until a phase shift of 90-deg was measured between the facing pressure 
and the cavity rear wall pressure. After resonance was determined pressure- 
phase data were obtained at that frequency and then at frequencies both above 
and below resonance. Data were obtained over the range of frequencies at sound 
pressure levels of 135, 140, 150, anci 160 db. 

The above procedure was repeated for each series of tests, but before 
each series the aperture diameter was increased by remachining the original 
aperture. Before the fifth series of tests, the individual resonator was 
transformed into a quarter-wave tube by increasing the aperture diameter to 
the cavity diameter and setting the backing depth equal 4. 64 in. , the theoretical 
quarter-wave length for a frequency of 730 Hz. In table II is a summary of the 
resonator dimensions for each of the experiments. 
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Table II. 

Dimensions of Acoustic Devices Used 
Verification Tests 

in First Series of 


Aperture 

Backing 

Theoretical 

Test 

Diameter, 

Cavity Depth, 

Resonant Frequency. 

Series 

in. 

in. 

Hz 

1 

0.75 

0.78 

750 

2 

1.0 

1. 16 

790 

3 

1. 5 

2. 22 

675 

4 

2.5 

3.53 

760 

5 

2.9 

4.64 ' 

730 


The acoustic resistance and reactance data were analyzed by computing 
the resistance coefficient, K, from 


K = R Cf/p u (55) 

and the effective aperture length coefficient, 6, from 

6 = ( t - t)/d ' (56) 

and comparing the results to the existing theory. 

For the resonators, as expected, the resistance coefficient was, on the 
average, the same as that of arrays, i. e. , 0.37; however, the quarter-wave 
tube was found to have a K value twice that of the array. This is because in a 
quarter-wave tube that is subjected to high incident sound fields a jet can form 
only during that part of the cycle in which there is flow out of the cavity whereas 
in a resonator a jet can form during both parts of the cycle. At high sound 
levels most of the acoustic energy is dissipated in the jet; thus, the tube will 
have approximately one-half the energy dissipation and, because dissipation and 
resistance are inversely proportional, twice the resistance of a resonator of the 
same aperture size. 

The experimentally determined effective aperture lengths were used in 
equation (56) to compute 5; an average value of 0. 66 was obtained which is in 
agreement with theory, equation (14). The averged value of 6, and a K value 
of 0. 37, was then used with the theory of Section II to compute resonator 
coefficients for the resonators. Two types of off-resonant calculations were 
made; the first by assuming the resistance was a constant equal to the resonant 
value, and the other, as explained in Section II, by solving the impedance 
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equation using an iterative procedure. Typical results are shown in figure 8. 

In the vicinity of resonance, either type of calculation agrees with the experiment; 
however, at frequencies far from resonance the iterative procedure, as expected, 
produced better results. 

A K value of 0. 75 and the iterative procedure was used to compute resonator 
coefficients for comparison with the quarter-wave tube data. As shown in 
figure 9 the agreement between theory and experiment was good. 



Figure 8. Comparison of Experimental Data from DF 90638 

Single Resonator With Nonarray Design 
Theories 


Figure 9. Comparison of Experimental Data from 
Quarter-Wave Tube With Design Theory 


DF 90639 
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It was concluded that if the resistance coefficient and the effective length 
correction factor can be specified for a given device, the nonarray theory can 
be used to accurately predict the resonator coefficient at resonance. Futhermore, 
the off-resonance coefficients can be predicted with sufficient accuracy for design 
purposes by either assuming constant resistance or employing an iterative 
solution to the impedance equation. 

C. LIMITATION EXPERIMENTS 

The limitation experiments were conducted to investigate the dimensional 
limitation of acoustic absorbing devices. Included was an investigation of the 
transition from Helmholtz resonator to quarter-wave tube, acoustic cavity 
experiments, aperture diameter and length experiments and resonator inter- 
action experiments. 

1. Resonator to Quarter-Wave Tube Transition Experiments 

A series of cold flow te.sts was conducted to study the transition between 
single resonators and quarter-wave tubes. The reverberation chamber was used 
as a sound source and the pressure-phase impedance technique was used to obtain 
data. Test sections consisted of two different resonators with adjustable backing 
depths; one resonator had a cavity diameter of 1.38 in. and the other 2. 90 in. 

Data were taken at seven different frequencies and at resonance; at each frequency, 
pressure-phase data were recorded at sound pressure levels of 150, 157, and 
165 db. 

An analysis was conducted of the impedance data obtained from each 
element to determine the limits for application of the design theory as the aper- 
ture of a single resonator is progressively enlarged until it becomes a quarter- 
wave tube, and to confirm the results obtained from the quarter-wave tube used 
in the previous series of experiments. 

It was found that there is essentially no transition between the single 
resonator and the A/4 tube; as shown in figure 10 the resistance coefficient 
averaged 0. 37 even when the ratio of aperture-to-cavity diameter approaches 
90%. Apparently only a small discontinuity in the aperture of a quarter-wave 
tube is all that is necessary to cause jet formation during both parts of the cycle. 
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Figure, 10. Comparison of Resistance Data from DF 61303 

Single Resonators and Quarter-Wave 
Tubes With Results from Previous 
Array Experiments 

The effective length correction factor data were found to be in agreement 
with theory. In figure 11 the data are shown with a line representing the theory, 
equation (14); it is worthy of note that the quarter-wave tube correction factors 
can be predicted with satisfactory accuracy. 

2. Acoustic Cavity Experiments 

Acoustic cavities are herein defined as those devices that do not cause 
the incident sound wave to pass through aperture(s) or other restrictions, e. g. , 
porous facings or screens, before entering the cavity. The type of acoustic 
cavity used for instability suppression has a cavity depth that is near the 
quarter-wave length of the incident wave. 

<, To study the energy absorbing capabilities of acoustic cavities, cold flow 
tests of three quarter-wave tubes and three quarter-wave slots were conducted 
using the high frequency impedance apparatus (Reference 2). The quarter-wave 
tubes had internal diameters of 0. 243 in. , 0. 305 in. , and 0. 430 in. and a depth 
of 1. 0 in. , which corresponds to the quarter-wave length for a frequency of 
3390 Hz. The quarter-wave slots (figure 12) each had the same cross-sectional 
areas as the quarter-wave tubes, and also had a cavity depth of 1. 0 in. The 
tests were conducted with each cavity by taking pressure-phase data at sound 
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levels from 140 to 160 db at resonance and also taking data over a range of 
frequencies from 1800 to 4000 Hz at a constant sound level of 150 db. 



Figure 11. Comparison of Effective Length DF 90640 

Correction Factor Data With 
Theory 






Figure 12. Quarter-Wave Slots Used in Acoustic FD 61304 

Cavity Experiments 
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All of the acoustic cavities were found to have fundamental frequencies 
lower than that corresponding to their quarter wave lengths. The ratio of the 
measured fundamental frequency to theoretical was approximately 0. 92 for 
the slots, and slightly lower for the tubes (~0. 9). Quarter-wave cavities 
have frequencies of resonance that are somewhat lower than theory because 
the air column in the device oscillates and thus has an overall or effective 
length that is slightly longer than the cavity depth. For the cavities used in 
these experiments, the effective length could be approximately determined 

by 

L = 1. 1L 
e 

where 

L = the actual cavity depth 

However, as shown previously, the effective length correction factor can be 
predicted using theory; therefore, a more accurate method of predicting the 
fundamental frequency of acoustic cavities would be to use equation (14) to 
compute 6 (for acoustic cavities, t = 0). Then use 


L = L + 6d 
e 

and 

f - c 
o 4L e 


(57) 

(58) 


to compute the frequency of resonance. If the device is not circular, an 
alternative form of equation (57) is 

-L = L + 0. 96 v/T (1-0.71 v/cT). (59) 

6 C 

Comparisons of the acoustic power absorbed by the slots and tubes showed 
that at resonance the energy absorbed by slots was greater than that absorbed 
by tubes even though each had identical open areas and backing depths. (See 
figure 13. ) In addition, it was found that the slots had a lower resistance 
coefficient (average value of K was 0. 53) than tubes for a given sound level, 
and the slots had average particle velocities approximately twice those of tubes. 

From further analysis it was concluded that the difference in peak power 
absorption between tubes and slots is caused by differences in the local particle 
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velocities at the edges of the two devices. Because the tube aperture is 
symmetrical, it has local particle velocities that are uniform around the perim- 
eter. In the slot the smaller dimension controls the local acoustic flow through 
the aperture, and because this dimension is always less than the tube diameter 
causes the acoustic flow to be more constricted than in the tube resulting in 
higher local particle velocities around the slot periphery for a particular 
incident sound pressure level. Energy absorption is a strong function of the 
local velocities, hence, the net result of the constriction is to cause the slot 
to have higher peak absolution. 



Figure 13. Comparison of the Power Absorbed by a DF 90641 
Quarter-Wave Tube and a Quarter-Wave 
Slot 


The bandwidth or off-resonance absorbing characteristics of the two devices 
were also different; as shown in figure 13, absorption with the slot changed with 
frequency more rapidly than with the tube indicating that the Q value for the slot 
was higher. These results are to be expected since the Q values of most 
acoustic devices are inversely proportional to resistance. 

The analysis of data from the acoustic cavity experiments was concluded 
by comparing the acoustic power absorbed by the three sets of quarter-wave 
tubes and slots with the theoretical absorption of single resonators having the 
same aperture areas and cavity volumes; results for one set are shown in 
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figure 14. In every case it was found that for the same sound level and frequency 
of resonance the equivalent single Helmholtz -type resonator will absorb more 
acoustic power than either a slot or a tube. This is because, as explained 
previously, the resonator dissipates energy in the two jets that form during 
each cycle; but with quarter-wave devices a jet can form only during that part 
of the cycle when there is flow out of the cavity, therefore less energy can be 
dissipated. 



Figure 14. Comparison of Acoustic Absorption Data DF 90642 
from a Quarter-Wave Tube and a Quarter- 
Wave Slot With the Theoretical Absorption 
of a Single Resonator 

From the results of this analysis, it is concluded that acoustic cavities 
are less desirable than resonators for use as rocket combustion stabilizing 
devices because: (1) for the same cavity volume and aperture area more energy 
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can be dissipated by a resonator; (2) the inherent narrow band absorbing char- 
acteristics and the difficulty in predicting the fundamental frequency make the 
tuning of acoustic cavities much more critical than that required for resonators. 
Even if a cavity is precisely tuned for the instability in a given motor and purged 
with gas of known properties to provide constant sonic velocity, little or no 
absorption will be obtained if the instability occurs at a slightly different 
frequency (which often happens). 

3. Aperture Diameter Experiments 

The objective of this experiment was to determine the diameter limitation 
of a resonator aperture. A series of pressure-phase experiments was conducted 
using the reverberation chamber as a sound source. A single resonator with an 
aperture diameter of 1. 25 in. and an adjustable cavity depth of 2. 9 in. diameter 
was used. The cavity depth was decreased in 0. 2 in. increments starting at 
1.2 in. ; data were recorded at the corresponding frequencies of resonance. As 
the backing depth was decreased the resonance frequency increased causing the 
wavelength to approach the aperture diameter. 

The acoustic data were first analyzed using the ratio of the aperture 
diameter to wavelength. The resistance coefficient was calculated from 
equation (55) and the effective length coefficient from equation (56); results are 
shown in figure 15. The resistance coefficient was constant until the aperture 
ratio exceeded 0. 3, then increased rapidly while the reactance, as represented 
by the effective length coefficient, increased slightly. 

It is known from previous work. Reference 2, that for operation in the 
nonlinear regime the particle velocity must be greater than 60 ft/sec. Examina- 
tion of the diameter limitation data revealed that the majority of the orifice 
velocity values were less than 60 ft/sec; thus a second series of diameter limi- 
tation experiments was conducted to obtain data at higher particle velocities. 

The tests were made with facing samples of 0.375 and 0. 182 in. diameter and 
backing cavity depths varying from 0. 90 to 0. 05 in. The data were again 
recorded at the resonance frequency of each configuration as the backing depth 
was decreased. The results are shown in figure 16; the data from the earlier 
experiment with orifice particle velocities greater than 60 ft/sec are also in- 
cluded. An increase in resistance coefficient was again noted for the values of 
4 d/A greater than 0. 3; however, the results also indicated that the resistance 
coefficients for the 0. 182-in. diameter experiments began to increase at 
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4 d/x = 0. 15. A correlation was then made (figure 17), between the resistance 
coefficient and the backing cavity depth to aperture diameter ratio. The 
resistance coefficient was found to be constant until the ratio of backing depth 
to aperture diameter decreased below 0. 8. A second correlation was made 
to combine the effect of aperture diameter, wavelength, and cavity depth on the 
resistance coefficient. The results are presented in figure 18, where the resis- 
tance coefficient as a function of 4 d/X L is shown. The data indicate that for 
values of 4 d/x L in excess of 1.0 the resistance coefficient deviates significantly 
from the value used in the design theory (0. 37). 

A third series of experiments was conducted using the high frequency 
impedance tube to verify the aperture diameter limitation. The test was conducted 
using a single aperture having a 0.213 in. diameter and a thickness of 0.125 in. 
The diameter was chosen so that wall interaction with the source side of the 
impedance tube would not affect the results. In addition, to prevent the high- 
frequency fundamental modes of the cavity from introducing error, the tests 
were conducted without a cavity. The results are shown in figure 19; results 
of the previous experiments are also included. Data for 4 d/X in excess of 0.4 
could not be obtained because of frequency and sound pressure limitations of 
the impedance tube. 



Figure 15. Comparison of the Resistance and Effective DF 90643 
Length Coefficients for the First Series of 
Aperture Diameter Limitation Tests 
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Figure 16. Comparison of the Resistance and Effective DF 90644 
Length Coefficient as a Function of Diameter 
to Wavelength Ratio 



• LO t.o 3.0 4.0 S.0 

BACKING DEPTH TO DIAMETER RATIO, L/d 


Figure 17. Correlation of Resistance Coefficient Data DF 90645 
With Backing Depth to Aperture Diameter 
Ratio 
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18. Correlation of Resistance Coefficient Data DF 90646 
With Diameter, Wavelength, and Cavity 
Depth 



Figure 19. Comparison of Resistance Coefficient Data DF 90647 
for Apertures With and Without Cavity as a 
Function of the Diameter Wavelength Data 
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The results indicated considerable resistance coefficient scatter was^ 
present for values of the diameter to wavelength ratio in excess of 0. 3 for the 

— s \ 

experiment without a cavity; however, the average value of the coefficient 
between 4 d/A =0.3 and 0. 4 was 0. 55, which is consistent with data from the 
earlier experiments. 

From the results of the limitation experiments it is recommended that 
apertures of single resonators be designed so that the ratio of aperture diameter 
to the product of backing depth and quarter-wave length is less than 1.0 at 
resonance (all dimensions in inches). 

4. Aperture Length Experiments 

The objective of the experiments was to determine the length limitations 
of resonator apertures. The high frequency impedance tube was used as a 
sound source and both the standing wave and the pres sure -phase impedance 
techniques was used to obtain data. Two series of experiments were conducted: 

(1) data were obtained at resonance for a configuration with apertures shorter 
than or equal to the quarter-wavelength, and (2) data were obtained with an 
aperture longer than quarter wavelength. 

a. Apertures Shorter Than or Equal to the Quarter-Wavelength 

An analysis was conducted of the impedance data obtained to determine 
the limits for application of the design theory as the aperture length was increased. 

It was found that the resistance coefficient did not change significantly from the 
predicted value of 0.37 as the thickness was increased from 0.5 to 3.35 in. 

(See figure 20. ) However, the frequency of resonance predicted by 


f o 27 r Y L/ < 

was found to be increasingly in error as the aperture length was increased. 

The error should be expected since equation (60) is obtained from the 
reactance equation at resonance 

X=pw/.ff £L = 0 i 

c u> L 

that was simplified from 


= pc tan(~j -pc<7 c cot 
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by the restriction that WL/c and 0)1/ c be small, so that cot o>L/c = c/coL and 
tan (j)£/c = o)H/ c . As shown in figure 21 this assumption and thus the simplified 
equation is accurate for values of o)i/c less than 0.4; however, for larger values, 
which is the case for the longer aperture configurations tested, the error becomes 
significant. 



Figure 20. Effect of Thickness on the Resistance DF 90648 

Coefficient 
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A comparison of the resonance frequency from the impedance data and the 
calculated data using both equation ( 61 ) and ( 62 ) are shown in table in. 


Table III. Results at Resonance for Apertures Smaller Than or 
Equal to Quarter-Wavelength 


Aperture 

Length, 

in. 

Resonance Frequency 
Impedance Data, 

Hz 

Resonance Frequency, 
Equation 61, 

Hz 

Resonance Frequency 
Equation 62, 

Hz 

0.5 

911 

957 

915 

1.5 

923 

1009 

936 

2.0 

909 

1088 

891 

2.5 

912 

1232 

909 

3.0 

910 

1590 

913 

3. 15 

913 

1740 

910 

3.35 

710 

1062 

733 

3.35 

704 

1012 

716 


It was concluded from the analysis that the design theory of Section II can be used 
with reasonable accuracy for values of jji/c and ooL/c less than or equal to 0.4, 
for other values equation (62) should be used in the calculation of the frequency of 
resonance and the acoustic reactance. 

b. Aperture Longer Than a Quarter-Wave Length 

An analysis was conducted of the impedance data to determine the effects 
of apertures longer than a quarter-wave length on absorbing characteristics. 

The data were obtained for a resonator having a 3. 25 in. long aperture with a 
backing cavity depth of 0. 262 in. , and an aperture diameter of 0. 18 in. The 
data were obtained at frequencies from 400 to 6000 Hz. Since the values of 
u)L/c and u)l/c were at times greater than 0. 4 for the frequency range of 
interest, equation (62) was used to calculate the absorption characteristics of 
the resonator, i.e., solve for the resonant frequency, fq, by letting X = 0 and 
use the equation to compute the reactance at other frequencies. As expected, 
it was found that in addition to the predicted fundamental mode resonance also 
occurred at several higher frequencies. The resistance coefficient was found 
to be approximately 0. 37 for the fundamental and for all harmonics. Figures 22, 

23, and 24 show comparisons of the experimental data with the theory for reactance, 
absorption coefficient, and power absorbed, each as a function of frequency. 
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Figure 22. Comparison of Theoretical and Experimental DF 90649 
Reactance for the Fundamental Mode and 
Harmonics of a Single Resonator 



Figure 23. Comparison of Theoretical and Experimental DF 90650 
Absorption Coefficients for the Fundamental 
Mode and the Harmonics of a Single Resonator 
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Figure 24. Comparison of Theoretical and Experimental DF 90651 
Power Absorbed for the Fundamental Mode 
and the Harmonics of a Single Resonator 

It was concluded that apertures longer than a quarter-wave length could be 
used in the design of acoustic resonators only if equation (62) is used to calculate 
the acoustic characteristics of the device. This procedure was followed in the 
formulation of the design computer program of Appendix C. The program will 
predict both fundamental and harmonic absorbing characteristics if any occur 
within the specified range of frequencies. 

5. Resonator Interaction Experiments 

Designers of acoustical absorbers may, for a given aperture area, use a 
single aperture or two (or more) smaller apertures. Equations (14) and (45), 
show in theory, that if smaller apertures are used the effective length will also 
be smaller, resulting in a higher frequency of resonance; thus, to maintain the 
same frequency of resonance the cavity volume must be increased. 

A series of experiments was conducted to investigate the effect of aperture 
interaction on resonance using a given total aperture area and various numbers of 
apertures all connected to a central, unpartitioned cavity. The high frequency 
impedance apparatus was used to perform the experiments; both pressure-phase 
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and standing wave data were obtained. Five resonator facings having 1, 2, 4, 6, 

2 

and 25 apertures were used. Each facing had a total open area of 0. 026 in. , a 

3 

cavity volume of 0. 028 in. , and facing thickness of 0. 1 in. The single aperture 
was tested first; the absorption was measured at resonance. The absorption of 
the remaining four resonators was then measured at two frequencies, one cor- 
responding to the resonant frequency of the single resonator and at the resonant 
frequency of the particular resonator. Results are listed in table IV. 


Table IV. Results From Single and Multiple 
Aperture Experiments 


Number of 
Apertures 

Resonance 
Frequency, Hz 

■Acoustic Power Absorbed, ft-lbf/sec 

at Resonance 

at 1580 Hz 

1 

1580 

0.52 

0. 52 

2 

1620 

0.52 

0.52 

4 

1780 

0.52 

0.52 

6 

1810 

0. 52 

0.50 

25 

1940 

0.52 

0.42 


It was found that within the limits of experimental error, all the resonators 
except the 25 aperture facing had the same absorption at resonance, and as the 
number of apertures was increased, the frequency of resonance increased. Thus, 
as theory predicts, to maintain resonance at the original frequency, i. e. , 1580 Hz, 
the cavity volume would have to be increased. These results show that if the reso- 
nator is not cavity volume limited, the number of apertures used to produce a 
given open area is immaterial. However, if the resonator is volume limited, 
highest absorption for a given cavity volume will be obtained with the fewest 
number of apertures. 

The 25 aperture facing had significantly lower abosrption at 1580 Hz than 
any other facing because of the usual decrease in absorption that occurs at 
frequencies far from resonance, and because of the small diameter of the aper- 
tures (0.037 in.). It is believed that as the aperture diameter is decreased 
below approximately 0. 050 in. the boundary layer begins to fill a significant per- 
centage of the open area. The average oscillatory velocity in the boundary layer 
is lower than the main stream particle velocity, thus less energy dissipation 
takes place and a decrease in the bandwidth characteristics occurs. For this 
reason, perhaps 0. 050 in. , is a realistic lower limit for aperture diameter. \ 

Experiments were next conducted to determine how a single resonator is 
affected when a second identical but completely separate resonator is brought 
into the near vicinity. Each resonator had a single aperture with a diameter of 
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0. 125 in. , a cavity diameter of 0. 305 in. , and a backing depth of 0. 875 in. The 
experiments were conducted using the high-frequency impedance apparatus; data 
were obtained from standing wave and phase measurements, i. e. , the resonant 
frequency was determined by varying the driver input frequency until a phase 
angle of 90 deg was obtained between the cavity and facing microphones. Standing 
wave data were then taken at resonance over a range of sound levels from 145 to 
170 db. 

It was found, as shown in table V, that introduction of the second element 
had essentially no effect on the frequency of resonance. 


Table V. Results of Interaction Experiments 


Sound Pressure Frequency of Resonance, Hz 

at Facing, db One Resonator Two Resonators 


144. 6 

1941 

1940 

149.8 

1944 

1944 

154.9 

1954 

1954 

159.4 

1961 

1961 

164. 9 

1974 

.1977 

169.8 

— 

2009 


The theory predicts that because neither the cavity open area ratio nor the 
aperture diameter were changed in the experiment, the frequency of resonance 
should have remained constant. The slight change in f Q with sound pressure, 
while certainly within the accuracy demonstrated by the correlation of figure 11, 
is not predictable by theory. For this reason it is suggested that nonarray ab- 
sorbers, like liner arrays, be designed for good off-resonant absorption charac- 
teristics. 

Comparison of other data from the experiment showed excellent agreement. 
The power absorbed data are shown in figure 25 and the particle velocity data are 
shown in figure 26. 

j 

The final series of interaction experiments was conducted to determine if 
placing the aperture of an absorbing device in a corner, e. g. , in a rocket 
chamber at the intersection of the chamber wall and the injector face, would 
affect its acoustical characteristics. Two resonators having four apertures of 
the same diameter and length connected to a common, unpartitioned cavity were 
fabricated for testing in the high-frequency impedance tube. The first, the con- 
trol resonator, had the apertures arranged so that each was as far removed from 
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the sidewall of the tube and from the other apertures as possible; the apertures 
of the second were located as close to the sidewalls as possible (figure 27). 
Each resonator was tested over a range of frequencies from 1550 to 1900 Hz at 
a sound level of 160 db. 



Figure 25. Results of Interaction Experiments With DF 90652 
Identical Resonators 
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Figure 26. Comparison of Aperture Particle Velocity DF 90653 
from Interaction Experiments With 
Theory 



Figure 27. Types of Facings Used in Corner Inter- FD 62348 
action Experiment 
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Two effects of placing the apertures in the corners were noted: (1) a 2. 5% 
increase in resonance from 1785 to 1830 Hz, and (2) more significantly, an in- 
crease in the power absorbed (figure 28). The small change in resonance is 
believed to be within the limits of experimental accuracy, and therefore, cannot 
be attributed to the presence of the corner walls, but, the increase in absorption 
can. Analysis of the absorption data revealed that the increase was produced 
primarily by a slight increase in particle velocity. Apparently, the corner walls 
produced a funneling effect, which resulted in slightly higher acoustic velocities 
during both halves of the acoustic cycle. 



Figure 28. Absorption Data from Corner Interaction DF 90654 
Experiment 

For the design analysis of resonators that are to be placed within one 
diameter of a wall or other structure, e. g. , injector face baffle, it is suggested 
that a value of 0. 34 be used for the resistance coefficient, which will cause the 
predicted absorption and particle velocity to be increased by the corresponding 
amount observed in the above experiments. 


■) 
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SECTION IV 

APPLICATION OF THEORY -SAMPLE CALCULATIONS 


The following calculations illustrate a method for determining the optimum 
physical dimensions of a nonarray absorbing device using the theory of Section II 
and the cold flow experimental results of Section III. The device will be single 
resonator; as shown in Section III, these devices are more desirable than any 
other type of nonarray absorber because of their inherently greater energy 
dissipation and broad-band absorbing characteristics. It is assumed that the 
resonator will be located in the wall of a rocket chamber so that the aperture 
is near the injector faceplane, the most effective location for such a device, 
and that the combustor will be operating at the following conditions: 

Chamber Pressure, P = 60 psia 

Combustion Gas Temperature, T c = 6000°R 

Frequency of Oscillations to be Damped, f 0 = 6000 Hz 

Gas Constant of Combustion Products, T = 45.7 ft lbf/lb m °R 

Ratio Specific Heats, 7=1.2 


Past experience (Reference 5) has shown that the aperture gas temperatures 
are significantly lower than the combustion gas temperature even if no liner 
cooling or purge flows are used. Although it is extremely difficult to predict 
the aperture gas temperatures in an untested motor, a reasonable assumption 
based on past experience is that the aperture gas temperature is one-half the 
combustion gas temperature, i. e. , 


T = T = 3000 °R 

3. Z C 


Using perfect gas relations the density 


and sonic velocity in the aperture 


j 


are 



(60) (144) 
(45. 7) (3000) 


= 0. 063 lb m /ft 3 


c = V 7 grT a = \/ 1. 2(32. 2)(45. 7)(3000) = 2300 ft/sec 


A. RESONANCE CALCULATIONS 

The fundamental frequency of the resonator will be the same as the 
frequency of the oscillations to be damped. At that frequency the wavelength 
is 

X Q = c/f 0 = 2300/6000 = 0.383 ft. 
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the wave number is 

k = 27T f / c = T 51 = 16. 41 ft" 1 
° A o 

and the critical spacing is 

b = 0. 4 \ Q = 1. 84 in. 

For maximum effectiveness, as many resonators as possible should be installed 
in the chamber; the above spacing is the minimum distance between the centers 
of adjacent apertures with separate resonator volumes. If for any reason the 
spacing must be closer than the minimum, the array theory (Reference 2) must 
be used for the design. 

1. Aperture Area 

For the zero flow case, i.e. , negligible net velocity past and/or through 
the aperture, the optimum aperture area is found using equation (33). 



By definition B= KCf where K is the resistance coefficient (0. 37 for circular 
apertures*) and Cf is the aperture discharge coefficient. For this case it is 
assumed that the aperture will have slightly rounded edges so that Cf = 1; 
then 

B = 0.37(1)" 2 = 0.37. 

As suggested in Reference 5, the resonator is designed for the maximum 
pressure amplitude allowed for combustion to be considered stable, usually 5% 
of the static pressure. Using 5% for this example we have pj/P = 0. 05 and the 
optimum area is 



For a single circular aperture 

d . - , /4(0. 0029)11442 . „ 73 ln , 

opt y ir V T 


*For other values of K see page 2. 
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If the optimum diameter is too large for the particular application it is possible 
to use two or more smaller apertures with a common cavity to obtain the required 
open area because, as shown in Section in, the number of apertures used to 
produce a resonator of given total open area with a common unpartioned cavity 
is immaterial unless the device is cavity volume- limited. If cavity volume is 
limited, highest absorption will be obtained with the fewest number of apertures. 

If multiple apertures are used it is recommended that all be contained within a 
surface area no larger than x Q /27r. 

The above restraints will cause, in some cases, applications in which it 
will not be possible to have aperture areas as large as required for optimum 
absorption. Then, a smaller area should be used and the resonator coefficient 
evaluated using the following procedure. For this example assume that a single 
aperture of diameter 0, 5 in. can be installed. The aperture area is 

2 

A = 7rd 2 /4 =■ = 0- 00136 ft 2 

and the dimensionless radiation resistance, computed using equation (30), is 


gk 2 A/2T =(16.41)V 00136) = 0 „583 
r 2T 

Compute the dimensionless particle velocity, which for zero flow is also the 
dimensionless dissipation resistance, using equation (29) 



= 0.1492 

The ratio of the two resistances, f} Q is 

0 = R/R =m/e = 0.149/0.0576 = 2.56 

o r o r 

and the resonator coefficient is computed using equation (11), which for resonance 
(X = 0) becomes 
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If the optimum aperture area had been used in the above calculations, 0 Q = 1 
and e = 1 would have resulted. Note that a 53% reduction in the aperture area 
decreased the resonator coefficient by only 19%. For the remainder of the 
computations necesjsary to design the resonator the optimum aperture area will 
be used, i. e. , d = 0. 73 in. 

2. Cavity Backing Depth 

The cavity backing depth is determined using the results of the Section III 
Limitation Experiments. From those experiments it was found that 


4d ... . -1 

< 1 m. 


*o L 


is a limiting condition, which as in this case, if the diameter and wavelength 
are known, can be used to determine the minimum backing depth, L m - n ; e. g. , 
rearranging the above equation produces 


< L 


hence. 


J rmn 


4d(l in. ) 


Substituting the known values (d and x q ) gives 

t 4(0.73) _ 

mm (0. 3 83) (12) 


r 


3. Aperture Length and Cavity Open Area Ratio 


The last two physical dimensions of the resonator may be obtained by 
first assuming a value for the aperture effective length, l, and then using the 
reactance equation (62) at resonance to determine the cavity open area ratio, 
<7 C . The aperture length, t, can then be computed from the effective length 
equation that applies for the particular device being designed. 

At resonance, equation (62) can be written as 

/2Tf 0 L\ 

*“ 1 — ) 
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Assuming a value of 0.35 in. for l and using L = L m _ n = 0. 635 produces 

= [ 27T(6000)(0.35)1 [27 r(6000)(0. 635) 

c |_ (2300)(12) J 2300(12) 

= tan (0.478) tan (0. 87) 

= (0. 518)(1. 188) = 0. 613 

It must be remembered that <J C is the area ratio between the aperture cross 
section and that of the cavity, not between the aperture and the chamber surface 
area. If the cavity has a circular cross section the diameter will be 

(l c ' V^e r 

J 

= 0. 73/v/O. 616 
= 0. 93 in. 

For this case (a resonator with zero flow) the applicable effective length 
equation is l = t + 6 d where 6 is computed from equation (14) 

5 = 0. 85(1 - 0.71ja~) 

1 / 

= 0. 85(1 - 0. 71 v/0. 616) = 0. 377 

Then, by rearranging the effective length equation, £= t + 5d, to obtain _ 

t = i- <5d 

= 0. 35 - 0. 377(0. 73) = 0. 075 in. 

The effective length equations for use with flow or for quarter-wave cavities 
are found on pages 9 and 28, respectively. 

B. OFF-RESONANCE CALCULATIONS 

The absorbing characteristics of the resonator must now be determined 
by computing the resonator coefficient at frequencies other than resonance. 

v o 

The quality factor is first computed using equation (40) 

© = x °A 

A(3 0 + l) 

_ 0.383(0,35) 

0. 0029(2)(12) 

= 1. 926 
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Assuming the resistance constant, the resonator coefficient can then be calcu- 
lated using equation (42) 

4 ^o 1 

(0 o + i) 2 i + Q 2 (r-i/n 2 

where f is the frequency ratio f/f Q - 

In table VI are the resonator coefficients calculated using this equation for 
frequencies from 5000 to 7000 Hz and in figure 29 is a comparison of the calcu- 
lated resonator coefficients with those computed for the example case using the 
computer program of Appendix C. The only difference between the two sets of 
calculations is the computer program does not use the constant resistance as- 
sumption to compute the off-resonance absorbing characteristics, (See page 11. ) 


Table VI. Calculated Values of Resonator Coefficient 
Based on Constant Resistance Assumption 


Frequency, Hz 

Frequency 
Ratio, f/f Q 

Resonator 

Coefficient 

5000 

0. 833 

0.67 

5200 

0.867 

0.77 

5400 

0.90 

0.86 

5600 

0.933 

0.93 

; 5800 

0.967 

0.98 

6000 

1.00 

1.00 

6200 

1. 03 

0.98 

6400 

1. 07 

0. 94 

6600 

1. 10 

0.88 

6800 

1. 13 

0. 81 

7000 

1. 167 

0.74 


C. FLOW EFFECTS 


A flow of gas past or through the apertures can, if the Mach number is high 
enough, cause an increase in the dissipation resistance and a change in the res- 
onator coefficient. 

1. Past-Flow Resonance' Calculations 

The following calculations show a method of evaluating the flow effects 
assuming a net flow past the aperture with a Mach number, (M ) of 0. 15. 
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Figure 29. Comparison of Computed Resonator FD 62364 

Coefficients for Example Case 

The Mach number of the oscillatory flow in the aperture at resonance is 
found using the definitions 



where m Q is computed using equation (29). Again using the optimized aperture 
area in equation (30) 



kA opt /27r 

(16. 4) 2 (0.0029)/27r) 


= 0.124 
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which when substituted in equation (29) gives 



0. 124 . i llO. 124 V 2 . 2(0. 37) (0. 05) 

" " 2 V\ 2 ) 1.2 

= 0. 124 

and an aperture Mach number of 


U o _ m o _ 0.124 _ 

— - or - °- 335 

To determine if the past-flow can be ignored the left side of inequality (45) 
is evaluated: 

2 

A 1C /aAAa\ ^ 

8 


M p / T c \ = 0. 15 / 6QQo \ 2 = , 

uj/c \T a ) 0.335 l 3000 1 


Because, as explained in Section II, the resulting value is greater than 1.5, the 
effects of flow cannot be ignored. The dimensionless dissipation resistance 
with flow is computed using equation (46). 

2 


R 

_R 

pc 


BM 


ft) •♦(*) 

0.37(0.15) / 6000\ 2 . 0.37 /n 
= J + _ (0.335) 


= 0.136 


The value of 0 is 
P 

R 

0 = 

P R r 


(Rp/Po) 


0. 136 
0. 124 


= 1.1 


which gives a resonator coefficient of 


6 = 
P - 



= 4(1- . 1) 

(2. I) 2 


0.997 


Thus, in this case, the flow has almost no effect on the resonator coefficient. 
If the coefficient had been reduced below approximately 0. 5 a new optimum 
aperture area should be found using equation (48) and the design procedure of 
subsection A repeated. 
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Because a different effective length equation is used with past-flow the 
aperture length must be recomputed. The effective length with past-flow, 
equation (44), can be written as /p = t + 0. 375 5d where all variables, 
except t, are the same as the zero flow case; therefore, 

t = l - 0.375 5 d 
P 

= 0.35 - 0.375(0.377) (0.73) 

= 0. 245 in. 

Hence, the most significant effect of the past-flow is the increase in aperture 
thickness (from 0. 075 to 0. 245 in. ) required to cause resonance to occur at the 
same frequency. 

2. Through-Flow, Resonance Calculations^ 

If there is a net flow of Mach number M^, through the aperture such that 
M t < u q /2c, the effect on absorbing characteristics can be ignored. If not, a 
procedure similar to that necessary with past-flow is followed except that 
equations (50) and (51) are used to compute the resistance and optimum area, 
respectively, and the effective length equation for zero flow, equation (14), is 
used. 

3. Off-Resonance Calculations 

The off-resonance absorbing characteristics with either type of flow are 
identical to those of the zero-flow case with one exception: with past-flow the 
quality factor, equation (40), is computed using the effective length from 
equation (44). 
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APPENDIX A 

REVERBERATION TEST TECHNIQUES 


A reverberation room is a highly sound- reflective room in which special 
care has been taken to make the sound field as diffuse as possible. Reverbera- 
tion tests are used to measure the absorption characteristics of acoustical mate-r 
rials. As the samples are relatively large and the sound is incident at random 
angles, the measurement of absorption is made, in many cases, under conditions 
closer to those of an actual installation than those in tests for impedance and 
absorption of acoustical materials by the impedance tube method. 

The sound that persists and decays in an enclosed space after the sound 
source has stopped is called reverberation. The rate of decay, D, in decibels 
per second is given by (Reference 4) 


D = 


1. 087 a c 
V 


(A-l) , 


where a is defined as the total absorptive power of the room of volume V ex- 
pressed as 



(A-2) 


with a ^ as the respective absorption coefficient representing the fraction of 
randomly incident energy absorbed by different surfaces of area, Sj, forming 
the interior walls of the reverberation chamber. 


In accordance with the original proposal of Sabine (Reference 9), who 
defined the reverberation time, Tgp, as the time required for the level of the 
sound in the room to decay to 60 db, then 



55.2 V 
c 

a c 


(A-3) 


In acoustical engineering it is customary to express Ve in cubic ft, 

and a in Sabines, where the Sabine is an absorption unit representing a surface 

O 

capable of absorbing sound at the same rate as does 1 ft of perfectly absorbing 
surface. Introducing for air, c - 1125 ft/sec, equation (A-3) becomes the 
Sabine reverberation equation 


0. 049 V 
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According to ASTM Standards (Reference 10), the Sabine equation is based on 
certain assumptions that must be fulfilled if measurements are to have meaning. 
These assumptions are: the sound field is diffuse before and during decay; the 
decay in decibels is essentially linear; and the measured decay rate is independent 
of the position of microphone or test sample. 

If an acoustical material having an absorptive surface, S, is suspended 
inside the reverberation chamber, a new reverberation time TgQ will equal 
(Reference 4) 



0.049 V 
c 

a + Sa 

e 


(A-5) 


where a Q is the effective absorption coefficient of the material being investi- 
gated. If equations (A-4) and (A-5) are combined to eliminate a, then 


°- 049 v c (_i_ 
e “ S \60 “ T 60 


(A-6) 


In recent years (Reference 11), there has been a tendency to standardize 
the quantity measured in reverberation chambers as the decay rate of sound 
pressure level in decibels per second rather than the reverberation time. Under 
such circumstances, it is advantageous to use equation (A-3) to eliminate Tg 0 
and TgQ from equation (A-6) and, thereby, obtain 


a 


e 


0.00081 V 

^ (D» - D) 


(A-7) 


where D' is the measured decay rate when the material being investigated is 
present in the reverberation chamber, and D is that measured when it is absent. 

For some acoustical elements, e.g. , quarter-wave tubes, slots, and non- 
array resonators, the surface area can only be arbitrarily defined, hence, a 
measure of total absorption is more meaningful. The total absorption, W, of 
the element can be obtained from reverberation data using the equations 

W = 0.049 V (A-8) 

C \ oO oO/ 


and 


W = 0.00081 V c (D' - D) 


(A-9) 
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The principal limitations of the reverberation room method are found at 
low frequencies (below 250 Hz), where the rooms are too small, and at fre- 
quencies above 4000 Hz, where absorption by the test media itself limits the 
accuracy. The works of Knudsen (Reference 12) on absorption caused by water 
vapor in the air pointed out the reason for the short reverberation times at high 
frequencies. Since water vapor absorption may be different when the empty 
room reverberation tests are made than they are when absorbing room tests 
are made, it can be a serious cause of inaccuracy, in addition to errors inherent 
in measurements of rapid decays. Consequently, steps should be taken to pro- 
vide a constant humidity for tests made in reverberation chambers. In addition, 
diffraction or area corrections are important for small porous test samples 
having high absorption. Parkinson (Reference 13) discussed experimental 
evidence of the diffraction problem, and, in 1939, Chrisler (Reference 14) gave 
the basis of empirical area correction used by the National Bureau of Standards. 

Comparisons of reverberation and calculated absorption coefficients have 
been made (Reference 15), in which the calculated values were based on impedance 
tube data. Close agreement is possible (figure A-l) for absorbers designed so that 
air can flow only normal to the surface. These absorbers are of the resonant-type, 
consisting of perforated panels with air space backing and are identical to full- 
array liners used in rocket chambers. Although excellent agreement exists be- 
tween the reverberation and impedance tube technique of measuring absorption 
characteristics, the reverberation method has a definite advantage in that non- 
array devices of any shape or configuration can easily be tested. 

Reverberation chambers, according to ASTM Standards (Reference 10), 
should be constructed of massive, nonresonant materials. The average absorp- 
tion coefficient of the chamber walls, determined by empty room measurements, 
should be less than 0.05 at all frequencies of interest, after allowance has been 
made for atmospheric absorption. The room must be isolated sufficiently from 
outside noises and vibrations to keep the ambient noise level 15 db below the 
lowest level at which measurements are made. A pentagonal chamber configura- 
tion is recommended in MIL-STD-810B, Method 515. A chamber of this configura- 
tion is best proportioned to produce, as closely as possible, a diffuse sound field, 
with a sound energy density that is very nearly uniform throughout the enclosure 
as required for use of Sabine’s equation in calculating absorption. 
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100 1000 

FREQUENCY - Hz 


Figure A-l. Comparison of Absorption Coefficients FD 40091 
Obtained from Reverberation and Im- 
pedance Tube Tests 

Test specimen size should be limited to a maximum value not to exceed 
10% of the chamber volume. Excessively small specimens should also be 
avoided to minimize the effects of diffraction (References 12, 13, and 14). The 
test specimen may be placed on the floor for larger chambers (over 8000 ft 3 ), 
but for smaller chambers it is necessary to suspend the specimen from the 
ceiling to ensure that the sound is randomly incident on the surface of the speci- 
men. 

The sound source used to drive the chambers should be one or more high 
intensity speakers mounted preferably in the trihedral corners of rectangular 
rooms. It may be necessary in oddly shaped and small chambers to mount the 
speakers outside the chamber and to introduce the test signal through exponen- 
tial horns mounted through the side walls of the chamber. 

Various experimenters (References 16 and 17) have tried to determine the 
best test signals to use when performing reverberation experiments. Pure 
tones, frequency modulated warble tones, combinations of pure tones (multi- 
tones), bands of pink noise, and pistol-shot noise have all been investigated. 

Use of a pure tone signal is not recommended because of a number of normal 
modes of vibration will be excited and, during the decay period, each of these 
normal modes assumes its own decay rate, resulting in beats that can cause 
irregular overall decay rates. (See figure A-2. ) No general agreement has 


r 
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been reached as to which signal is best but, based on comments of Schultz 
(Reference 18), the use of either a warble tone or bands of white noise is suit- 
able for accurate measurements of reverberation time. 




Figure A-2. Decay Curves With Double Slopes Pro- FD 40088 
duced by Normal Modes of Vibration 
With Different Decay Rates 

Typical electronic equipment used in reverberation measurements is shown 
in a block diagram of figure A-3. A tracking generator produces a pure tone 
signal, which is then input into a frequency modulator to produce a warble tone 
output. The warble tone signal is then amplified and introduced into the rever- 
beration chamber through two high- intensity loudspeakers. Four microphones 
are located throughout the chamber to determine if the sound field is diffuse and 
to provide decay rates for various locations in the chamber as recommended by 
ASTM Standards (Reference 10). The microphone signal is then passed through 
a narrow band filter center at the same frequency as the pure tone signal. A 
logarithmic potentiometer is then used to convert the filter microphone signal to 
a dc output, which is recorded on a graphic recorder in the form of decibels 
versus time. A high speed oscillograph recorder could be used to give an 
accurate measurement of a 20-db drop occurring in a matter of tenths of a 
second. The graphic recorder thus, provides a paper-tape record of the decay 
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of sound in the reverberation room (figure A-4). Repeated measurements of 
the decay rate under the same conditions should be expected to show some 
scatter. In making reverberation tests, therefore, it is necessary to make 
several decay curves for each mean frequency tested and also for each position 
of the microphone. From the data an average reverberation decay rate is ob- 
tained and used in equation (A-7) or equation (A-9) to calculate the absorption 
characteristics of the test specimen. 



Figure A-3. Schematic of Instrumentation Used for FD 40084 
Measurement of Decay Time in a Re- 
verberation Test 



Figure A-4. Typical Graphic Record of the Decay FD 40090A 
of Sound in a Reverberation Chamber 
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APPENDIX B 

GENERAL PRESSURE-PHASE DATA REDUCTION EQUATIONS 


The following equations were used to determine the impedance of acoustical 

devices from measurements of the sound pressure directly in front of the aperture, 

Pj_, in the cavity, P 2 , and phase difference, 4 > , between P-^ and P 2 . The device 

must consist of a cavity of uniform cross sectional area, A , and arbitrary length, 

c 

L, exposed to an incident sound wave through an aperture of cross sectional area A 
no larger than that of the cavity. 


Impedance 



(sin^ + i cos^) 


Resistance (real part of Z) 


(B— 1) 


g c P C 

sin k L 



sin4> 


Reactance (imaginary part of Z) 


X = 


~ a c pC 
sin k L 


p- J COS0 
2 > 


Particle Velocity in Aperture (peak) 
s/T g P 2 sin k L 
°c<> c " - 


(B-2) 


(B-3) 


(B-4) 


Power Absorbed 

A P P„ g (sin k L) sin<*> 

C 1 z 


w = 


2 pc 


Resistance coefficient defined by R = K pu C^' 
D / u W \ 2 

K = 


-2 


g c cC I 
sin k L 


sin<*> 


v/i~ g 

Effective Aperture Length 
<7 


l = k Ti nkL (cos k L " < P 1 /P 2> C0S< P> 


(B-5) 


(B-6) 


(B-7) 
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APPENDIX C 

DESIGN COMPUTER PROGRAM 

In this Appendix is contained a digital computer program that computes the 
theoretical acoustic absorbing characteristics of nonarray devices. The program 
was written in FORTRAN IV for the IBM 1130 computer. The contents appear in 
the following order: 

1. Input Format 

2. Output Symbol Nomenclature 

3. Sample' Input 

4. Source Program Listing 

5. Output for Sample Cases 

The sample input sheet may be used as input for a program test case; 
results identical to those of the Example Cases should be obtained. 
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. 1. Input Format 


CARD 

COLS 

VARIABLE 

DESCRIPTION 

I 

01-05 

NCASE 

NO. OF COMPLETE SETS OF INPUT 

2 

01-80 

TITLE 

ANY TITLE INFORMATION 

3 

01-10 

PRESS 

STATIC PRESSURE (PSIA) 


11-20 

GASCN 

APERTURE GAS CONSTANT ( LBF-FT/LBM-R ) 


21-30 

TEMP 

APERTURE GAS TEMPERATURE (RANKINE) 


31-40 

TEMPC 

COMBUSTION GAS TEMPERATURE (RANKINE) 


41-50 

GAMA 

RATIO SPECIFIC HEATS. APERTURE GAS 


51-60 

CF 

APERTURE DISCHARGE COEFFICIENT 


61-70 

T 

APERTURE LENGTH (INCH) 


NOTE- 

- FOR QUARTER WAVE TUBE USE INPUT T-0.0 


71-80 

L 

BACKING CAVITY DEPTH (INCH) 

4 

01-10 

D 

APERTURE DIAMETER (INCH) 


11-20 

SIGC 

APERTURE-TO-CAVITY OPEN AREA RATIO 


21-30 

FM IN 

MINIMUM FREQUENCY OF INTEREST (HZ) 


31-40 

DELF 

FREQUENCY INCREMENT (HZ) 


41-50 

FMAX 

MAXIMUM FREQUENCY (HZ) 


51-60 

SPLN 

MIN PRESSURE RATIO. PEAK INCIDENT AMPL 1 TUDE-TO-STAT 1 1 


61-70 

DELP 

PRESSURE RATIO INCREMENT 


71-80 

SPLM 

MAX. PRESSURE RATIO 

5 

01-10 

MP 

MACH NO. OF FLOW PAST APERTURE 


11-20 

MT 

MACH NO. OF FLOW THRU APERTURE 


21-30 

FO 

FREQUENCY OF RESONANCE (HZ)-IF BLANK THE FUNDAMENTAL 


FREQ OF THE ABSORBER WILL BE COMPUTED. 


NOTE ALL CARDS (EXCEPT 1) MUST BE INPUT NCASE TIMES 


2. Output Symbol Nomenclature 

VARIABLE DESCRIPTION 


FREQ 

RFREQ 

RESIS 

REACT 

RADRX 

UORFX 

BETAX 

RCOEF 

POWAB 

CRITA 


FREQUENCY (HZ) 

RATIO OF FREQUENCY TO RESONANT FREQ. 

DISSIPATION RESISTANCE/ (DENSITY * SONIC VELOCITY 
REACTANCE /(DENSITY * SONIC VELOCITY) 

RADIATION RESISTANCE /(DENSITY * SONIC VELOCITY) 
ORIFICE PARTICLE VELOCITY (FT/SEC) 

RESISTANCE RATIO. D I SS I PAT I ON-TO-RAD I AT ION 
RESINATOR COEFFICIENT 
POWER ABSORBED (FT-LBS/SEC) 

RATIO. ORIFICE DIAMETER TO QUARTER WAVE- 
LENGTH AND BACKING CAVITY DEPTH ll./IN.) 
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4. Source Program Listing 


SUBROUTINE SOLVM < FUNCT * YACT . X . TOL .DELA. YCAL ) 
DEL-DELA 
UPAS ■ 1 
1*1 

10 YCAL«FUNCT< X) 

IF ( ABS ( (YCAL-YACT) /YCAL >-TOL>l. 1*2 

2 GO TO (3.4) »KPAS 

3 GO TO (5.6) .1 

5 1=2 

ystor-ycal-yact 

XSTOR=X 

IF (YCAL-YACT) 7.1.8 

7 KFLAG-1 

15 X=X+OEL 
GO TO 10 

8 KFLAG-2 

16 X=X-0EL 
GO TO 10 

6 KPAS= 2 
Y=YCAL-YACT 
IF(Y*YSTOR 130.1. 11 

11 I F ( ABS ( YSTOR l-ABS ( Y) ) 12.1.13 

12 DEL = -DEL 

13 YSTOR = Y 
XSTOR = X 

GO TO <15.161 .KFLAG 

4 Y= YCAL-YACT 

IF( Y*YSTOR) 17.1 .18 

18 YSTOR ■ Y 
XSTOR-X 

GO TO (19.20) .KFLAG 

19 X=X+DEL 
GO TO 10 

20 X=X— DEL 
GO TO 10 

17 DEL-DEL/10. 

30 GO TO (21.22) .XFLAG 

21 X-XSTOR+DEL 
GO TO. 10 

22 X-XSTOR-DEL 
GO TO 10 

1 CONTINUE 
RETURN 
END 

// DUP 

•STORE W5 UA SOLVM 

// JOB T 
// FOR 

•ONE WORD INTEGERS 

FUNCTION FUNKL(X) 

REAL L 

COMMON L.XLEFF.SIG 
Z=L/12. 

W=XLEFF/12. 

FUNKL=SIG*COS( X*Z )*COS( X*W) /SIN (X*2 ) /SIN ( X*W) 
RETURN 
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END 

// OUP 

♦STORE WS UA FUNKL 

// FOR 

♦ONE WORD INTEGERS 

FUNCTION FUNKUIX) 

REAL L 

COMMON L.XLEFF.S1G.RRD.XD.PIPP 

FUNKU ■ X*X* ( ( X+RRD ) **2 + XD**2) / PIPP**2 

RETURN 

END 

// DUP 

♦STORE WS UA FUNKU 

// FOR 

♦IOCStCARD. 1132 PRINTER) 

♦ONE WORD INTEGERS 
** DECK 13172 

EXTERNAL FUNKL , FUNKU 

REAL L.KP.MP.MT.KO.MO.KL.KOL.KT 

DIMENSION TITLEI20) 

COMMON L.XLEFF.S1G.RADRX.REACT tPIPP 
99 FORMAT (15) 

100 FORMAT(BEIO.O) 

READ (2*99) NCASE 
DO 406 MWH= 1 tNCASE 
READ (2*101) TITLE 

101 FORMAT ( 20A4 ) 

READ ( 2 » 100 ) PRESS. GASCN.TEMP ♦ TEMPC • GAMA »CF . T . L 

READ ( 2 • 100 ) D.SIG . FM I N .DELF . FMAX . SPLN .DELP . SPLM 

READ ( 2 • 100 ) MP.MT.FO 

P1«SPLN 

G=32« 1725 

PIO.1416 

PI2-6.2B32 

CF2« CF*CF 

XK2-.37 

YACT-1. 

TOL-.OOl 

DEN=144.0*PRESS/GASCN/TEMP 
SON I C = SQRT ( GAMA*GASCN*TEMP*G ) 

APARE°PI*(D**2) /4.0 
XLEFF=T+.85*D* 1 1. - . 71*SQRT l S I G ) ) 

I F ( FO ) 9 . 8 . 9 

9 ANGLE = PI2*FO*L/12« /SON I C 

XLEFF = SIG*SONIC/P I 2/F0*C0S( ANGLE) /SIN( ANGLE > * 1 2 • 
I F ( T ) 32.32.53 
3? XK2»2.*XK2 
GO TO 53 

8 IF(MP) 10.11.10 

10 XLEFF = T + .375*( .S5*D*( l.-.7*SQRT(SIG) > > 

11 IF (T) 5.1.5 

1 FO » S0NIC*12./4./(L+XLEFF> 

XK2*2 •* XK2 
GO TO 53 

5 FO • SORT ( SIG/L/XLEFF) *12.* S0NIC/PI2 
KP=PJ2*F0/S0NIC 
X=KP 

DELA=X/4.0 
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CALL SOLVM ( FUNKL * YACT , X , TOL . DELA . YCAL ) 

FO=X*SONIC/P ! 2 
53 KO= PI2«FO/SONIC 
303 WR I TE ( 3 » 20 1 ) TITLE 

201 FORMAT ( ' 1 ' 1X,20A4,//32X'GAS PROPERTIES'/) 

WR I TE < 3 « 202 ) PRESS • TEMP «GASCN. GAMA , DEN .SON IC 

202 FORMAT (lOX'STATlC PRESSURE ( PS I A ) ' FAB. 3/ 

110X' APERTURE GAS TEMP (DEG. R)'F44.1/ 

2 10X ' GAS CONSTANT < L8F*FT/LBM/R ) *F43. 3/ 

3 10X ' R AT 10 OF SPECIFIC HEATS'F47.3/ 

4 10X 1 GAS DENSITY < LBM/FT**3 ) ' F47.4/ 

5 10X • GAS SONIC VELOC I TY ( FT/SEC )' F44. 1 ) 

WR I TE ( 3 • 203 ) TEMPC.MP.MT 

203 FORMAT ( 10X ' COMBUST 1 ON GAS TEMP (DEG R) ' F43.1 / 
llOX'MACH NO OF FLOW PAST APERTURE' F41.3 / 

2 10X ' MACH NO OF FLOW THRU APERTURE* F41.3 // ) 

WRITE(3*205) SIG.CF.T.D.L 

205 FORMAT ( ' '31X'LlNER PARAMETERS'// 

UOX'CAVITY OPEN AREA RATI0'F48.4/ 

2 10X ' FLOW C0EFFIC!ENT'F54.3/ 

310X* APERTURE LENGTH < IN) ' F50.3/ 

410X' APERTURE DIAMETER <IN)'F48.3/ 

5 10X ' CAV I TY BACKING DEPTH (IN)'F45.3> 

WR I TE ( 3 * 206 ) XLEFF • AP ARE . XK2 

206 FORMATIIOX'EFFECTIVE APERTURE LENGTH ( IN) 'F40.3/ 

1 10X ' APERTURE AREA ( I N**2 ) ' F49.4/ 

210X'RESISTANCE COEFFICIENT' F48.3 // ) 

RRFO « K0*K0*APARE/144./PI2 '~ 

B= XX2/CF2 

PIPP »2.* B*P1 /GAMA 

MO ■ SORT ( RRF0*RRF0/4.+P I PP I -RRFO/2. 

UORFO » MO*SON I C/B 
RESFO - MO 
MO ■ MO/B 

IF (MT-MO/2.) 22.21.21 

21 RESFO » MT/CF**3 
GO TO 23 

22 I F ( MP*50N IC/U0RF0*(TEMPC/TEMP)**2 -1.5) 23.23.24 
24 RESFO « B*MP/3.*1 TEMPC/TEMP)**2 -t-M0/2.*B 

23 BETAO ■ RESFO/RRFO 

EPSLO «4.*BETA0 / ( BETAO+1. )**2 
QUALO«SONIC*XLEFF*12.0/APARE/(BETAO+1.0) /FO 
PI ■ P1*PRESS/S0RTI2. > 

PID «20.*.4343*AL0G(PI/29.E-10) 

PIP«PI*SQRT(2. ) 

ENAB0«,5*DEN*APARE/144.*( SON I C*MO ) **3*B/G 
WRITE (3.207) P I O.P IP .UORFU .RESFO, RRFO 

207 FORMAT ( • '20X'LINER ACOUSTIC PARAMETERS AT RESONANCE*// 
HOX'RMS AMPLITUDE OF INCIDENT PRESSURE (DECIBELS)' F25.2 / 

2 10X ' PEAK AMPLITUDE OF INCIDENT PRESSURE ( PS 1 A ) ' F28.5 / 

3 10X * OR I F ICE VELOCITY AT RESONANCE ( FT /SF.C I ' F32 . 1 / 
410X'DISSIPATI0N RESISTANCE (DIMENSIONLESS)' F32.3 / 

51 OX 'RAD I AT ION RESISTANCE ( 0 IMENS I ONLESS ) ' F34.3 ) 

WRITE (3.208) BETAO .EPSLO .OUALO, ENABO , FO 
20fl FORMAT ( lOX'RESI STANCE RATIO. DI SS I PAT I ON-TO-RADI AT I ON ' F28 . 3/ 
1 10X • RESONATOR COEFFICIENT AT RESONANCE ' F36. 3/ 

210X ' QUALITY FACTOR ' F56. 3/ 

3 10X ' ENERGY ABSORBED ( FT*LBF/ SEC ) ' E42 .4/ 
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4 10X ' RESONAN T FREQUENT (H2)'F48.1 ) 

FREQ-FMIN 
7 WRITE (3*6) 

6 FORMAT ( 1H1 I 

WRITE <3*2001 TITLE 
200 FORMAT ( 1X.20A4.//) 

LINE = 1 
WRITE (3*209) 

209 FORMAT <11X' FREQ '4X'RFREQ'4X' REACT '6X'RESIS'4X' RADRX' 3X' UORFX' 3X 
l'BETAX'3X'RC0EF'4X' POWA8 ' 6X ' CR I TA ' / ) 

603 CONTINUE 

RFREQ-FREQ/FO 

CR I TA-D-FREQ/L /SON IC/3 « 

RADRX « RRFO*RFREQ**2 
KOL-KO-L/ 1 2 « 

KL-KOL*RFREQ 

KT-KO*RFREQ*XL£FF/12. 

REACT«<SlG*COSIKL)/SIN(KL)-SIN(KT)/COS<KT> ) 

X=MO*B 

DELA-X/4.0 

CALL SOLVM (FUNKU.YACT.X.TOL.DELA.YCAD 
UORFX “ X*SONIC/B 
RESIS « X 

IF (MT-UORFX/SONIC/2. ) 25*26*26 

26 RESIS - MT/CF**3 
GO TO 27 

25 TR- (TEMPC / TEMP) **2 

I F ( MP*SONIC/UORFX*TR- 1.5) 27.27.28 
28 RESIS- B*MP/3.*TR+ X/2. 

27 BETAX - RESIS / RADRX 

RCOEF - 4.0*BETAX/( (BETAX ♦ 1. 01**2 + ( REACT/RADRX ) **2 > 
ENABX«£nAB0*RESIS/RESF0*(U0RFX/U0RF0>**2 

WRITE (3*211 ) FREQ. RFREO. REACT. RES1 S. RADRX. UORFX, BET AX. RCOEF *ENABX 
l.CRITA ^ 

LINE - LINE + 1 

211 FORMAT! F16.1. F8.3.2F10.4.F9.4, F8 . 2 . 2F8 . 3 . 2E12 . 4 ) 

I F ( FREQ-FMAX ) 601.602.602 

601 FREQ- FREQ+D ELF 

IF (LINE-50) 603.7,7 

602 IF (Pi - SPLM) 401. 406. 406 
401 P1-P1+DELP 

GO TO 303 
406 CONTINUE 
CALL EXIT 
END 

// XEQ 
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5. Output For Sample Cases 


example case* no flow 

GAS PROPERTIES 

STATIC PRESSURE (PSIA) 

APERTURE GAS TEMP (DEG. R) 

GAS CONSTANT ( LBF*FT/L8M/R ) 

RATIO OF SPECIFIC HEATS 
GAS DENSITY (LBM/FT**3) 

GAS SONIC VELOC ITY(FT/SEC) 

COMBUSTION GAS TEMP (DEG R) 

MACH NO OF FLOW PAST APERTURE 
MACH NO OF FLOW THRU APERTURE 


LINER PARAMETERS 

CAVITY OPEN AREA RATIO 
FLOW COEFFICIENT 
APERTURE LENGTH (IN) 

APERTURE DIAMETER ( IN) 

CAVITY BACKING DEPTH (IN) 

EFFECTIVE APERTURE LENGTH (IN) 

APERTURE AREA (IN**2) 

RESISTANCE COEFFICIENT 


LINER ACOUSTIC PARAMETERS AT RESONANCE 


RMS AMPLITUDE OF INCIDENT PRESSURE (DECIBELS) 177. 2B 
PEAK AMPLITUDE OF INCIDENT PRESSURE (PSIA) 2.99999 
ORIFICE VELOCITY AT RESONANCE (FT/SEC) 771.9 
DISSIPATION RESISTANCE (DIMENSIONLESS) 0.124 
RADIATION RESISTANCE (DIMENSIONLESS) 0.124 
RESISTANCE RATIO. D I SS IPAT ION-TO-RAD I AT ION 0.999 
RESONATOR COEFFICIENT AT RESONANCE 1.000 
QUALITY FACTOR 1.927 
ENERGY ABSORBED ( FT*LBF/SEC ) 0.4844E 03 
RESONANT FREQUENY (HZ) 6000.4 


0.6131 

1.000 

0.075 

0.730 

0.636 

0.350 

0.4185 

0.370 


60.000 

3000.0 
45.700 

1.200 

0.0630 

2300.6 

6000.0 
0.000 
0.000 


68 



Pratt& Whitney Aircraft 

PWA FR-4993 


EXAMPLE CASE * NO FLOW 


FREO 

RFREO 

REACT 

RESIS 

RADRX 

UORFX 

RETAX 

RCOEF 

POWAB 


CR1TA 


5000.0 

0.833 

0.2723 

0.0943 

0.0862 

586.65 

1.093 

0.304 

0.2126E 

03 

0.8315E 

00 

5100.0 

0.849 

0.2429 

0. 1000 

0.0897 

621.97 

1.114 

0.377 

0.2534E 

03 

0 • 848 IE 

00 

5200.0 

0,866 

0.2141 

0.1055 

0.0932 

656.13 

1.131 

0.461 

0.2975E 

03 

0.8647E 

00 

5300.0 

0.883 

0.1859 

0.1106 

0.0969 

688.16 

1.141 

0.552 

0.3432E 

03 

0.8813E 

00 

5400.0 

0.899 

0.1581 

0.1152 

0.1006 

716.72 

1.145 

0.647 

0.3877E 

03 

0.8980E 

00 

5500.0 

0.916 

0.1308 

0.1191 

0.1043 

740.65 

1.141 

0.741 

0.4279E 

03 

0.9146E 

00 

5600.0 

0.933 

0.1039 

0.1220 

0.10B2 

759.18 

1.128 

0.827 

0.4608E 

03 

0.9312E 

00 

5700.0 

0.949 

0.0773 

0. 1240 

0.1120 

771.53 

1.106 

0.900 

0.4837E 

03 

0.9479E 

00 

58C0.0 

0.966 

0.0512 

0.1250 

0.1160 

777.70 

1.077 

0.955 

0.4954E 

03 

0.9645E 

00 

5900.0 

0.983 

0.0254 

0. 1250 

0.1201 

777.70 

1.041 

0.988 

0.4954E 

03 

0 . 981 IE 

00 

6000.0 

0.999 

-C.00C0 

0.1241 

0.1242 

771.91 

0.999 

1.000 

0.4844E 

03 

O.9970E 

00 

6100.0 

1.016 

-0.0252 

0.1223 

0.1283 

760.91 

0.953 

0.989 

0.4640E 

03 

0. 1014E 

01 

6200.0 

1.033 

-C.Q502 

0,1197 

0.1326 

744.90 

0.903 

0.959 

0.4353E 

03 

0.1031E 

01 

6300.0 

1.049 

. *0.0749 

0.1166 

0.1369 

725.21 

o.esi 

0.913 

0.4017E 

03 

0.1047E 

01 

6400.0 

1.066 

-0.0994 

0.1129 

0.1413 

702*44 

0.799 

0.856 

0.3650E 

03 

0.1064E 

01 

6500,0 

1.083 

-0.1236 

0.1089 

0.1457 

677.35 

0.747 

0.792 

0.3273E 

03 

0.1080E 

01 

6600.0 

1.099 

-0.1477 

C. 1046 

0.1502 

650.72 

0.696 

0.724 

0.2902E 

03 

0.1097E 

01 

6700.0 

1.116 

•'-0.1716 

0. 1002 

0.1548 

623.32 

0.647 

0.656 

0.2550E 

03 

0.1114E 

01 

6800.0 

1.133 

-0.1954 

C. 0959 

0.1595 

596.30 

0.601 

0.591 

0.2233E 

03 

0.1130E 

01 

6900.0 

1.149 

-0.2190 

0.0915 

0.1642 

569.29 

0.557 

0.530 

0.1943E 

03 

0.1147E 

01 

7000.0 

1.166 

-0.2424 

0.0873 

0.1690 

543.43 

0.516 

0.474 

O'. 1690E 

03 

0.1164E 

01 
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example case* past flow 

GAS PROPERTIES 


STATIC PRESSURE (PSIA) 60.000 
APERTURE GAS TEMP (DEG. R) 3000.0 
GAS CONSTANT ( LBF*FT /LBM/R ) 45.700 
RATIO OF SPECIFIC HEATS 1.200 
GAS DENSITY <lBM/FT*»3> 0.0630 
GAS SONIC VELOC ITYIFT/SEC) 2300.6 
COMRUSTION GAS TEMP (DEG R> 6000.0 
MACH NO OF FLOW PAST APERTURE 0.150 
MACH NO OF FLOW THRU APERTURE 0.000 

LINER PARAMETERS 

CAVITY OPEN AREA RATIO ' 0.6131 
FLOW COEFFICIENT 1.000 
APERTURE LENGTH (IN) 0.245 
APERTURE DIAMETER (IN) 0.730 
CAVITY BACKING DEPTH (IN) 0.636 
EFFECTIVE APERTURE LENGTH (IN) 0.350 
APERTURE AREA (IN**.?) 0.4185 
RESISTANCE COEFFICIENT 0.370 

LINER ACOUSTIC PARAMETERS AT RESONANCE 

RMS AMPLITUDE OF INCIDENT PRESSURE (DECIBELS) 177.28 
PEAK AMPLITUDE OF INCIDENT PRESSURE (PSIA) 2.99999 
ORIFICE VELOCITY AT RESONANCE I FT/SEC) 771.6 
DISSIPATION RESISTANCE (DIMENSIONLESS) 0.136 
RADIATION RESISTANCE l DIMENSIONLESS) 0.124 
RESISTANCE RATIO. DISSIPATION-TO-RADIATION 1.093 
RESONATOR COEFFICIENT AT RESONANCE 0.997 
QUALITY FACTOR 1.837 
ENERGY ABSORBEO ( FT*LBF/SEC ) 0.4839E 03 
RESONANT FREGUENY (HZ) 6003.7 
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example case* past flow 


FREQ 

RFREQ 

REACT 

RES I S 

RADRX 

UORFX 

BETAX 

RCOEF 

POWAB 


CR1TA 


5000.0 

0.832 

0.2728 

0.1211 

0.0862 

586.24 

1.404 

0.355 

0.2487E 

03 

0.8315E 

00 

5100.0 

0.849 

0.2435 

0.1239 

0.0897 

621.16 

1.381 

0.423 

0.2856E 

03 

0.8481E 

00 

5200.0 

0.866 

0.2147 

0.1267 

0.0932 

655.50 

1.358 

0.500 

0.3252E 
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